Gauge invariant response functions in Algebraic Fermi liquids 
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A new method is developed that permits the simple evaluation of two-loop response functions for 
fermions coupled to a gauge field. We employ this method to study the gauge-invariant response 
functions in the Algebraic Fermi liquid, a non-Fermi liquid state proposed to describe the pseudogap 
phase in the QED3 theory of cuprate superconductors. The staggered spin susceptibility is found 
to exhibit a characteristic anomalous dimension exponent t]a, while other correlators show behavior 
consistent with the conservation laws imposed by the symmetries of the underlying theory. 



Low-energy effective theories of certain correlated elec- 
tronic systems are known to be equivalent to (2+1) di- 
mensional quantum electrodynamics (QED3) 0,0, [J 
0,0 . The latter describes N species of massless 'relativis- 
tic' Dirac fermions coupled to a massless gauge field a M . 
While the physics behind these different reincarnations 
of QED3 varies from case to case, these theories are all 
of considerable general interest for the following reason: 
In the so-called symmetric phase of QED3, long-range in- 
teractions mediated by the massless gauge field produce 
characteristic anomalous correlations between electrons 
which decay on long length- and time-scales as nontriv- 
ial power laws. This symmetric phase has been termed 
variously as the Algebraic Spin [3 or Algebraic Fermi 
liquid (AFL) and embodies a unique realization of a non- 
Fermi liquid state of electronic matter in 2 spatial dimen- 
sions. 

The power law correlations are encoded in the fermion 
propagator of the theory, G(x) = ('F(x)*(0)) ~ :zH 2+ ''), 
where r\ — (4/3ir 2 N)(3t; — 2) is the anomalous dimen- 
sion exponent and £ > is a gauge fixing parameter. 
The above fermionic propagator is gauge-dependent and 
therefore it cannot represent the behavior of the phys- 
ical electron in the underlying theory. Much effort has 
gone into constructing and calculating t he p roper gauge- 
invariant electron propagator 0, 1 10| but the sit- 
uation remains unclear. The most natural candidate, 

Q(x — x') = (exp(i f a^ds^)^(x)^(x')), suffers from se- 
vere ultraviolet divergences due to the line integral of the 
gauge field and is meaningless in the absence of a phys- 
ically motivated regularization scheme 0. Attempts to 
evaluate Q directly have yielded an unphysical negative 
anomalous dimension pi IM lio|| . 

In this Letter we search for anomalous power law cor- 
relations in various physical response functions of QED3 
which represent susceptibilities and conductivities of the 
spin and charge degrees of freedom in the underlying the- 
ory. By construction these quantities describe gauge in- 
variant physical responses of the system to an external 
probe and therefore no difficulties arise as to the inter- 
pretation of the results. We evaluate the requisite corre- 
lators to leading nontrivial order in the 1/7V expansion 
[Tl|. The principal technical hurdle in any such calcu- 
lation is the evaluation of the vertex correction which is 



essential to preserve gauge invariance [12|. In the typi- 
cal treatment 0, 0] this leads to lengthy algebra due to 
overlapping singularities within Feynman diagrams. Here 
we devise a new method for isolating the leading diver- 
gent behavior of such diagrams which reduces the task 
to computing the trace of small number of Dirac gamma 
matrices. We test this new method against known re- 
sults and derive new results with essentially no extra ef- 
fort. We formulate a simple general rule for determining 
whether or not a given correlator exhibits any anomalous 
dimension and discuss this result in terms of a theorem 
that prohibits the correlators of conserved operators from 
acquiring an anomalous dimension 0, 0] . 

In what follows we focus on the QED3 theory of 
cuprate pseudogap as formulated in Refs. [(J |9( but our 
technique remains applicable to any other reincarnation 
of QED 3 H H H. In the former QED 3 emerges as 
a low-energy effective theory for the nodal topological 
fermions ^i(x), in a phase-disordered d-wave supercon- 
ductor (g?SC). The Lagrangian density is 
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where the gauge field a M encodes the topological frustra- 
tion encountered by fermions as they propagate through 
the "soup" of fluctuating unbound vortex-antivortex 
pairs, x = (r, r) denotes the space-time coordinate and 
7/i = 0,1,2) are 4x4 gamma matrices satisfying 
{7m> lv] = 25u!/. The bare dynamics of the gauge field is 
Maxwellian [(J 03- 113 with stiffness K^. 

In the above theory ^i(x) is a four component 
spinor describing the fermionic excitations at the Z-th pair 
of antipodal nodes of the underlying dSC. Its individual 
components are related to the original electron operators 
c„ through the singular gauge transformation 

c CT (r,r)=e^( r ^^(r,r), (2) 

detailed in Refs. 0, ||- The purpose of this transfor- 
mation is to "unwind" the phase ip(r, r) = (fii(r, r) + 
ifil(r,T) of the fluctuating SC order parameter A(r,r) = 
Aoe l¥, ( r,r ) in favor of coupling to the gauge field a p which 
is related to coarse-grained gradients of phases ip a . The 
difficulty in computing the gauge- invariant electron prop- 
agator Q{x — x') stems from the necessity to evaluate 
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FIG. 1: Feynman diagrams for susceptibility Eq. © to lead- 
ing order in 1/iV. Solid lines represent the bare fermion prop- 
agator, wavy lines represent the screened gauge-field propa- 
gator, and the solid squares stand for vertex T;. 



der in 1/N expansion we require 4 diagrams depicted 
in Fig. ^a-d). The wavy line represents the gauge field 
propagator that becomes universal at long wavelengths 
due to the screening by topological fermions [(|, 
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in the sense that Kp drops out and only reappears as an 
upper cutoff of the theory, A ~ 1/K. 

The bare bubble Fig. ^a) reads 



Xo (q)= f TriGoipjTtGoip-q)^ 

J V 
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the averages of such phase factors. The latter are dom- 
inated by short length-scale physics that is not properly 
described by the effective theory 

Here we focus on the charge and spin density correla- 
tors which do not suffer from the above problem. Con- 
sider in particular the spin susceptibility, 



dTe™ r (T T S*(T)S*(0)), (3) 



with (3 = 1/T and = / d 2 re lq r cr4(r) C(T (r) the 
^-component of the electron spin operator at wavevector 
q. The charge susceptibility \ c is defined as in Eq. © 
but with «Sq replaced by the charge density operator p q = 

J <i 2 re lqr J^a c o-( r ) c o-( r )- Both spin and charge densities 
are gauge singlets with respect to a^, i.e., they do not 
pick up any phase factors under the transformation 
Following Ref. Q we adopt the representation for the 
gamma matrices given by 7 M = a 3 ® (02, 01, — 03) with 
cr, the Pauli matrices. We further define 73 = <j\ <S> 1 and 
75 = 02 <8 1 and note that these anticommute with all 
7 M 's. With these definitions we may express the density 
operators as 



S z (q) 1 _ 
P{q) 



#i(p)rj*i(g-p), 



(4) 



where Go(p) =i>/p 2 = Pul^/p 2 is the free Dirac propaga- 
tor, and can be evaluated by standard methods (see e.g. 
Appendix B of Ref. |{|). One obtains 



Xo(q) = -Tr^r^r,]— 
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We observe that the structure of the bare bubble suscep- 
tibility is determined entirely by the commutation rela- 
tion of the vertex with the 7 matrices. We shall see that 
this property remains true for the higher order diagrams 
as well. 

Diagrams (b), (c) and (d) represent the leading 1/N 
corrections due to the gauge field and contain all the 
interesting physics. We note that their sum is gauge 
invariant. This follows since they represent correlators 
of gauge invariant density operators. One can explic- 
itly verify that this is so by employing the Ward identity 
Gq 1 (p + k) — Gq 1 (fc) = p M 7 M . We are thus free to fix the 
gauge and in what follows we adopt the Feynman gauge 
(£ = 1) in which D^(q) = (8/N)6fu,/q. Diagrams (b) 
and (c) contain the self-energy insertion and read 



X(b+c)(q) 



Tr[G (p)r i G (p - g)r ; G (p)S(p)] (9) 



where vertex Ti takes the form Tf = (170,270) an d 



= (— «7o72, *7o7i) for spin and charge densities with 
momentum transfer close to q = (0, 0), respectively, and 
Tf = (75,0), = (—7075,0) for momentum trans- 
fer near q = Qi m (tt,tt). We have switched to Eu- 
clidean notation with 3-momenta q — (qo, q), J denotes 

J d 3 p/(2ir) 3 and summation over I = 1,2 is assumed. 
The susceptibilities are then given as 

X(q) = f f <*i(p)r,¥ I (p + g)*i(|/)r,¥,(p' -<?)>. (5) 

We evaluate x(o) at T = by formally considering N/2 
identical copies of fields \&i and iff 2- Then to leading or- 



with S(p) = J fc D A(!/ (fc)7 M Go(p - k)j u =-Ti F ^>ln(A/p), 
and 



3ir 2 N 



(10) 



the fermion anomalous dimension in the Feynman gauge 
Inserting T*(p) into Eq. @ and combining with Eq. 
jSJ) one can express the divergent contribution as 

X(b+c)(q) ZVFXoiq) ln(A/g). (11) 

Diagram (d) contains the vertex correction and reads 
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X(d)(q) = JJtt([G (p - ?)r,G (p)]7^[Go(p - k^Gafr - k - q^D^k). (12) 

As mentioned above the direct evaluation of this type of integral presents a significant challenge [l3L fl4| . We are 
interested in the leading q — > behavior of x(q)- I n this limit the p integral has singularities as p — > and p — » fc, 
corresponding to divergences in the first and second term in the square brackets respectively. The main contribution 
to the integral therefore comes from the vicinity of these two points and we may evaluate it by expanding the regular 
part of the integrand at the singular point. Retaining only the leading term thus gives 

X{d) (q) ~ ^Tr([Gob-g)r ; G (p)h^Go(-A : )r ; Go(-fc- ( z)]7,)^W 

+ JJji^Goik - q^tGoik^iGoip - k^Goip - k - q)]^D^(k). (13) 



Performing a variable shift p — > p + k in the second term 
we may simplify the above expression as 

X(d)(«)=2Tb:(ni(g)jf[Go(p-g)r,Go(p)]) (14) 

with Q,(g) = f k7l ,{G (k)r l G {k-q)}'y 1/ D^(k). The last 
integral is again easy to evaluate and gives 



Mi(q) = [j^ v Ti^ v -/^]riF\D.(A/q). 



(15) 



We now notice that for vertices composed of products 
of gamma matrices (such as those entering the spin and 
charge densities defined above) it holds that 



7m 7^7*7/* = XT i, 



(16) 



where A is a number. In particular we shall encounter 
two types of vertices. Type-I vertices commute or an- 
ticommute with all 7 M 's (e.g. T t = 1,73,75,737s) and 
in this case A = 9. Type-II vertices anticommute with 
one or two of 7 M 's and commute with the rest (e.q. 
T( = 7o, 7o7i, ■ ■ • )• In this case A = 1. With this in- 
sight Eq. (|14f> becomes 



X ( d) (q) = 2rj F XHA/q) / r&[Go(p)r,Go(p-g)r,]. 
J P 

Combining this result with Eqs. (|11|) and Q we can write 
the result for the full susceptibility to 1/N order, 

X (q) = Xo(q)[l - 2vf(1 - A) ln(A/«)]. (17) 

This correction may be interpreted as the leading term 
of a power law [lll | , so that we have 

^ f A \ m 8(A-1) 
X(q)-X0(q)[-) , * = -^p (18) 

The anomalous dimension exponent r/ 4 is entirely de- 
termined by the algebraic properties of the vertex F; 



through Eq. ||TSJ|. In particular for type-I vertex (A = 9) 
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64 



(19) 



while for type-II vertex 774 = 0. 

Returning back to physics we see that the vertex for 
staggered spin susceptibility is type-I and will therefore 
exhibit nontrivial anomalous dimension exponent Eq. 
(JTDJ, in agreement with the result of Ref. Q obtained by 
laborious direct evaluation of the vertex correction. Uni- 
form spin susceptibility has type-II vertex and therefore 
does not acquire anomalous dimension from diagrams (b- 
d), again in agreement with |14| . In addition both charge 
susceptibilities are type-II and will not exhibit anomalous 
dimension. 

Finally we note that if coincides with one of the 
7 M 's then diagrams of the type shown in Fig. ^e) are 
nonvanishing and must be included in the leading order. 
One can show that resummation of diagrams (e) modifies 
the xo(q) ~ q behavior in Eq. (JSJ) to q 2 . This can be 
viewed as another type of anomalous dimension due to 
the coupling to the gauge field and we shall discuss it 
more fully elsewhere. 

The formalism we have developed allows us to study 
spin and charge conductivities in the AFL. These can 
be calculated through the Kubo formula as Cy (w) = 
-Imnjf(w)/w, where U\f(uj) = U i:j (iuj -> w + i5) with 



Uij(iu>) 



dre^(T Tjl (T)j 3 (0)) (20) 



the current-current correlation function. Indices i,j — 
1,2 label the spatial components of the 3-current j M . 
The spin current is given as j^(x) — ^ , ;(x)7 M 'I';(a;) while 
the electric current is given as = (*i7o72 v I'i — 

, F27o7i x I'2, 17172*1, -*27i72 v I / 2)- To leading 1/N 
order the computation of the correlator Q20(l involves the 
same diagrams depicted in Fig. ^ an d we can therefore 
simply adopt the result obtained for the susceptibility in 
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Eq. I|18[) . We see that the vertices involved in j s and 
are all type-II and therefore neither spin nor charge 
conductivity will exhibit anomalous dimension. 

We may conclude that of all quantities considered only 
the staggered spin susceptibility will bear the unique sig- 
nature of the AFL in that its frequency momentum and 
temperature dependence will be controlled by the anoma- 
lous dimension exponent 774 given in Eq. (|19fl . The other 
quantities will behave essentially as if the theory was un- 
correlated. 

The fact that some quantities remain essentially un- 
affected by strong long range correlations can be traced 
to a field theoretic theorem 0, 0] which states that 
the correlator of conserved currents (i.e. currents satis- 
fying dfj.jp, = 0) must exhibit a scaling dimension which 
agrees with its engineering dimension. Physically this 
means that some quantities are constrained by their con- 
servation laws to such a degree that correlations can- 
not alter their long-distance scaling behavior. As an 
example consider the spin current j*. It is a con- 
served current, <9 M j* = guaranteed by the gauge in- 
variance of the theory [i.e. invariance under local sym- 
metry ^/(a;) — > e ie ^ x '^i(x)]. As we have verified by ex- 
plicit calculation neither spin conductivity nor uniform 
spin susceptibility exhibit anomalous dimension (~ 1/N) 
beyond the xo ~ <Z 2 behavior mentioned above, in ac- 
cord with the theorem. Understanding the behavior of 
the charge conductivity is less straightforward as we do 
not expect the quasiparticle current to be conserved in 
a (phase-disordered) superconductor. Indeed there is no 
symmetry that would guarantee conservation. How- 
ever, the theory is known to posses a "crural" symmetry 
which produces two conserved cur- 
rents: jjp — ^;7 Al (j7375) v I'i (no sum over I). By virtue of 
the identity 707375 = —7172 these are related to the spa- 
tial components of j^, explaining the lack of anomalous 
dimension in the electrical conductivity. 

In the 4x4 representation of Dirac gamma matrices 
there are 4 type-I vertices listed below Eq. (|lt)|) . Thus 
there should be 3 other physical observables which ex- 
hibit anomalous scaling dimension. What are these? 



In turns out that 73 corresponds to the staggered spin 
susceptibility at the wave vector Qj = — Qi, while 1 
and 7375 correspond to susceptibilities to formation of 
subdominant (phase-incoherent) SC order in s and p 
channels, respectively. These quantities are members of 
the QED 3 chiral manifold [H G3, which is the mani- 
fold of broken-symmetry states occuring for N < N c w 
32/-7T 2 = 3.24 at zero temperature. This phenomenon 
is known as spontaneous chiral symmetry breaking [T^ |. 
The anomalous dimension exponent 774 that we find in 
the above susceptibilities appears to anticipate the tran- 
sition into the chiral symmetry broken phase. Indeed if 
we combine Xo(q) from Eq. (JSJ with Eq. I|18|) we find 
x(s) ~ ? 4 ; implying divergent q — > susceptibility for 
N < 64/3tt 2 = §JV C . 

In the context of the QED3 theory of cuprates [f| we 
expect the experimental manifestations of the nontrivial 
anomalous dimension to appear in the quantities inhab- 
iting the QED 3 chiral manifold, most prominently the 
staggered spin susceptibility that is directly measurable 
by neutron scattering. In addition, a similar nontrivial 
response should obtain in the charge channel near (0, 7f) 
as a consequence of enlarged chiral manifold discussed in 
[Isj . This particular response is unique to the present 
theory and will not be present in the SU(2) theory of 
Ref. Q. 

We conclude by observing that the physics of the Alge- 
braic Fermi liquid bears distinct similarity to that of ID 
Luttinger liquids. While the single particle properties are 
distinctly non-Fermi liquid like, only certain measurable 
physical responses exhibit the unique fingerprint of the 
AFL in the form of anomalous scaling dimension. The 
method developed here allows us to easily identify these 
quantities and helps in our search for situations where 
other quantities might develop anomalous behavior. 
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